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We study the extraction of the ground-state parameters from vacuum-to-vacuum correlators. 
We work in quantum-mechanical potential model which provides the only possibility to probe the 
reliability and the actual accuracy of this method: one obtains the bound-state parameters from 
the correlators by the standard procedures adopted in the method of sum rules and compares these 
results with the exact values calculated from the Schrodinger equation. We focus on the crucial 
ingredient of the method of sum rules — the effective continuum threshold — and propose a new 
algorithm to fix this quantity. In a quantum-mechanical model, our procedure leads to a remarkable 
improvement of the accuracy of the extracted ground-state parameters compared to the standard 
procedures adopted in the method and used in all previous applications of dispersive sum rules in 
QCD. The application of the proposed procedure in QCD promises a considerable increase of the 
accuracy of the extracted hadron parameters. 
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1. INTRODUCTION 

The extraction of the ground-state parameters from the operator product expansion (OPE) series for a relevant 
correlator is a cumbersome procedure: even if several terms of the OPE for the correlator are known precisely, the 
numerical procedures of the method of dispersive sum rules [l|, Q cannot determine the true exact value of the bound- 
state parameter. Instead, the method should provide the band of values such the true hadron parameter has a flat 
probability distribution within this band 1]. This band is a systematic, or intrinsic, sum-rule uncertainty. 

The method of sum rules in QCD contains a set of prescriptions (see e.g. [3, 4]) which are believed to provide such 
a systematic error. In QCD this, however, always remains a conjecture — it is impossible to prove that the range 
provided by the standard sum-rule procedures indeed contains the actual value of the bound-state parameter. 

The only possibility to acquire an unbiased judgement of the reliability of the error estimates in sum rules is to 
apply the method to a problem where the parameters of the theory may be fixed and the corresponding parameters of 
the ground state may be calculated independently and exactly. Presently, only quantum-mechanical potential models 
provide such a possibility. 

A simple harmonic-oscillator (HO) potential model [1, H, [1, 0] possesses the essential features of QCD — confinement 
and asymptotic freedom Q — and has the following advantages: (i) the bound-state parameters (masses, wave 
functions, form factors) are known precisely; (ii) direct analogues of the QCD correlators may be calculated exactly. 
(For a discussion of many aspects of sum rules in quantum mechanics we refer to [9I. [Tol fTTl fl3. fl3l [T^| 1 . 

Making use of this model, we have already studied the extraction of ground-state parameters from different types 
of correlators: namely, the ground-state decay constant from two-point vacuum-to- vacuum correlator Q|, the form 
factor from three-point vacuum-to- vacuum correlator [f| , and form factor from vacuum-to-hadron correlator § . We 
have demonstrated that the standard adopted procedures for obtaining the systematic errors do not work properly: 
for all types of correlators the true known value of the bound-state parameter was shown to lie outside the band 
obtained according to the standard criteria. These results give a solid ground to claim that also in QCD the actual 
accuracy of the method turns out to be much worse than the accuracy expected on the basis of the standard criteria 
applied. 

Thus, our results reported in [3, [H, @] mainly contained cautious messages concerning the application of sum rules 
to hadron properties. However, understanding the problem is the necessary first step in solving the problem. 

We have realized that the main origin of this failure of the method lies in an over-simplified model for hadron 
contnuum which is described as a perturbative contribution above a constant Borel-parameter independent effective 
continuum threshold. We have introduced the notion of the exact effective continuum threshold, which corresponds to 
the true bound-state parameters: in a HO model the true hadron parameters — decay constant and form factor — are 
known and the exact effective continuum thresholds for different correlators may be calculated. We have demonstrated 
that the exact effective continuum threhsold (i) is not a universal quantity and depends on the correlator considered 
(i.e. it is in general different for two-point and three-point vacuum-to- vacuum correlators), and (ii) depends on the 
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Borel parameter and, for the form-factor case, also on the momentum-transfer. 

Moreover, we have shown that the "Borel stability criterium" combined with the assumption of a Borel-parameter 
independent effective continuum threshold leads to the extraction of the wrong bound-state parameters. 

In view of these results a natural question arises: it is possible to formulate modified procedures which would 
improve the actual accuracy of hadron parameters extracted from sum rules, or the method is only suitable for 
obtaining numerical estimates with uncontrollable errors which may amount to 20-30%? 

We believe the answer to the first part of this question is positive. 

In a recent paper [tJ we proposed a new algorithm for extracting the parameters of the ground state. The simple 
idea of Q is to relax the standard assumption of a Borel-parameter independent effective continuum threshold, and 
allow for a Borel-parameter dependent quantity. In 0] we have illustrated the way the procedure works for the 
extraction of the form factor from three-point vacuum-to- vacuum correlator. 

In the present paper we develop this idea and provide details of its application to the extraction of hadron parameters 
from dispersive sum rules. We present for the first time the application of the new algorithm to the two-point correlator 
and present a detailed analysis of the extraction of the form factor from three-point correlator. The new algorithm 
is shown to lead to a significant improvement of the actual accuracy of the extracted bound-state parameters. Even 
though the full control over the systematic uncertainties, in the strict mathematical sense, is still not feasible (and 
presumably cannot be achieved in principle), our findings suggest that it may be possible to shed light on the actual 
accuracy of the method of sum-rules in QCD. 

The paper is organized as follows: Section [5] recalls the results for bound states in the HO model; Section [3] gives the 
OPE for the polarization operator. Section 0] presents detailed results for the vertex function and the corresponding 
OPE. Sum rules for the decay constant and the form factor are given in Section [5J Section [5] provides our numerical 
results. Finally, Section [7] presents our concluding remarks. 

2. HARMONIC-OSCILLATOR MODEL 

We consider a nonrelativistic model Hamiltonian H with a HO interaction potential V(r), r = \ f\: 

H = H + V(r), H =p 2 /2m, V(r) = mu 2 r 2 /2. (2.1) 

The full Green function G(E) = (H — E)~ x and the free Green function Gq(E) = (H — E)^ 1 are related by 

G-\E)-G^{E)=V. (2.2) 

The solution G(E) of this relation may be easily obtained as an expansion in powers of the interaction V: 

G(E) = G (E) - G (E)VG (E) + ■■■ . (2.3) 

In the HO model, all characteristics of the bound states are easily calculable. For instance, for the ground state (g), 
n = 0, one finds 

E K =-u, i? g ^ g (r = 0)| 2 = (— ) , F g (q) = exp(- g 2 /4mc.), (2.4) 

where the elastic form factor of the ground state is defined in terms of the ground-state wave function ip s according 
to 

F g (q) = (* s \J(q)\%) = J d S kiPl(k)iP s (k - q) = J d 3 r |^ g (f)|V^, q ee \q\, (2.5) 

and the current operator J(q) is given by the kernel 

(f'\J(q)\f) = eMiqr)5 {3) (r - r '). (2.6) 

3. POLARIZATION OPERATOR 



The quantum-mechanical analogue of the polarization operator has the form 

H(T) = (r) = 0| exp(-i/T)|r 4 = 0). 



(3.1) 
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This quantity is used in the sum-rule approach for the extraction of the wave function at the origin (i.e., of the decay 
constant) of the ground state [l[ . A detailed analysis of the corresponding procedure for the HO model can be found 
in For the HO potential, the analytic expression for n(T) is well-known 



\ 7T / 



nCO-fgT , ' w - (3-2) 

[2 sinn(wJ )J 



Apart from the overall factor, n(T) is a function of one parameter Tlo. 
The average energy of the polarization function is defined as 

E n (T) = -d T logU(T) = -u JC oth(u J T), d T = ^. (3.3) 

At T = both n(T) and En(T) diverge. For large values of T the contributions of the excited states to the correlator 
vanish and therefore En(T) — > E g for T — > oo. The deviation of the energy from i5 g at finite values of T measures 
the "contamination" of the correlator by the excited states. 

The ground-state contribution to the correlator has a simple form: 

n g (T) = (^) 3/2 exp(-^T). (3.4) 

The OPE series is the expansion of n(T) at small Euclidean time T: 

n-^^n^+j^ 1 +••■)■ <»» 

For n(T) this expansion is equivalent to the expansion in powers of the interaction to. 

The first term in (|3.5[) does not depend on the interaction and describes the free propagation of the constituents. 
Ho may be written as the spectral integral Q 

n (T) = J dze- zT p (z), p(z) = -^LVS, (3.6) 

with po(z) the known spectral density of the one- loop two-point Feynman diagram of the nonrelativistic field theory 
The rest of the series represents power corrections, which may be obtained just as the difference between the 
exact correlator and the free-propagation term: 

n powcr (T) = n(T)-n (T). (3.7) 

4. VERTEX FUNCTION 

The basic quantity for the extraction of the form factor in the method of dispersive sum rules is the correlator of 
three currents • The analogue of this quantity in quantum mechanics reads [l3| 

r(E 2 ,E 1 ,q) = {r f = 0\(H-E 2 )- 1 J(q)(H-E 1 )- 1 \r i =Q), q=\q\, (4.1) 

[with the operator J(q) defined in l|2.6p ] and its double Borel (Laplace) transform under E\ — ► t\ and E 2 — > r 2 

r{ T2 ,n,q) - (f f = 0\G(T 2 )J(q)G(n)\ri = 0), G(r) = exp(-iJr). (4.2) 

For large T\ and t 2 the correlator is dominated by the ground state: 

r(r 2 ,r 1 , 9 )^|^ g (r-0)| 2 e -^^+^F g ( (Z 2 ) + ... . (4.3) 

Let us notice the Ward identity which relates the vertex function at zero momentum to the polarization operator: 

r(T2 ) T 1 ,g = 0)=n(ri+r 2 ). (4.4) 

This expression follows directly from the current conservation relation 

J(q = 0) = 1. (4.5) 
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Fig. 1: Relative ground-state contribution to the energy E S (T, 
for several values of q. Solid (red) line: q — 0; dashed (green) 1 
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(b) 

?r(T,g) (a) and to the correlator T g {T,q) /T(T, q) (b) vs. T 
q — dash-dotted (blue) line: q — 



A. Exact T(T, q) in the HO model 

We obtain the exact T(t2, n, (?) by using the known expression for the Green function in configuration space 

(r ' = 0|G(r)|r) = H(r) exp tanh" 1 ^)) . (4.6) 

All necessary integrals are Gaussian and we easily derive an explicit expression for T{t2,ti, q). For our further 
investigation, we consider the vertex function for equal times T\ = T2 — ^T, which reads 

T(T,q) = II(T)exp f-^-tanh (^-)) . (4.7) 



Arauj \ 2 

The correlator T is a function of two dimensionless variables ojT and q 2 /mu>. Setting, without loss of generality, 
vn = lj we still have 2 dimensionless variables ojT and q/uj. This will be done later for the numerical analysis. 
The corresponding average energy is defined as follows 

E T (T,q) , -d T lo g T(T,q) = ^coth(.T) + ^ (1 + J^) - (4-8) 
The contribution of the ground state to the correlator, r g , in the HO model has the form 
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r«<r, 9 ) = (=) 3 '% XP (>)c Xp (-^) 



(4.9) 



The average energy Ey and the ground-state contribution T g vs. T for different values of q are shown in Fig. [T] 

First, notice that the on-set of the ground-state dominance in the correlator is shifted to later times with growing 
q. Second, if the ground-state energy E g is known (hereafter we will be discussing precisely this situation), the 
deviation of the average correlator energy from the ground-state energy may be used as an indicator of the ground- 
state contribution to the correlator. This indicator is, however, not very accurate: the relative contribution of the 
ground state to the correlator turns out to be systematically greater than the deviation of the average correlator 
energy from the ground-state energy. 

Figure [T] makes obvious the way in which the ground-state form factor may be extracted from the correlator T(T, q) 
known numerically (e.g., from the lattice): The correlator is dominated by the ground state at large values of T; so 
one may calculate the T- and g-dependent energy Er(T,q) (|4.8[) which exhibits a plateau at large T: E(T,q) — ► E g 
for any q. Making sure that one has already reached the plateau and that the correlator is saturated by the ground 
state, one obtains the form factor from the relation 

F s (q)= lim ± e ^ T r(T,g). (4.10) 

T— >oo Kg 
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B. OPEfor r(T,g) 

Let us now construct for T(T, q) the analogue of the OPE as used in the method of three-point sum rules in QCD. 
The corresponding procedure consists of two steps: First, we expand T in powers of u> 2 and obtain 

r(T j(Z ) = ^r 2n ( g ,rv 2 ". (4.n) 

n=0 

Explicitly, for the lowest terms one has 

,3/2 / q^T" 



= (sir) exp 



8 m 



The first term, To, corresponds to free propagation and does not depend on the confining potential. It may be written 
as the double spectral representation [15| 



r (T, q) = j dz 1 dz 2 e-^ lT e-^ 2T Ao{z l7 z 2 , q), A (zi, z 2 , q) = -£^8 ( ^ + z 2 - |^ - A Zl z 2 < j . (4.13) 



The Ward identity in nonrelativistic field theory relates to each other the three-point function at zero momentum 
transfer and the two-point function and leads to 

lim A (z, z', q) = S(z - z')p (z). (4.14) 

For the calculation of higher-order terms T 2n , n > 1, the explicit form of the confining potential is necessary. In 
the HO model it makes no problem to calculate all higher-order terms explicitly. Recall, however, that in QCD the 
situation is different: the confining potential is not known and therefore nonpcrturbative long-distance effects are 
parameterized as power corrections in terms of local condensates. Thus, for each quantity T 2n , n > 1, one has a 
power-series expansion in T. To keep the same track, we expand T 2n , n > 1, in powers of T. As the result of this 
procedure, the quantum-mechanical analogue of the OPE for T takes the form 



TopeC? 1 , q) = r (T, q) + r p0 wor(r, q) 



m \ 3 / 2 



r p ower(T\ q) — (2^7^) 



-w 2 T 2 + - T + uj 4 T 

4 24m I 480 1536m 2 



(4.15) 



We display here only terms up to 0(T 4 ) in r p0 wcr- For the practical extraction of the ground-state form factor we 
can retain any number of higher-order terms, which may be easily generated from the exact expression (|4.7p . 

It should be emphasized that the coefficients of each power of T n in the square brackets of (|4.15[) are polynomials 
in q 2 of order (n — 2). Therefore, if the momentum q increases, one needs to include more and more power corrections 
in order to achieve a certain accuracy of the truncated OPE series for T. In QCD, this implies the necessity to know 
and include condensates of higher and higher dimensions. Since only a few lowest-order condensates are known, in 
practice the applicability of three-point sum rules in QCD is restricted to the region of not too large q 2 and T. 

Figure [21 demonstrates the behaviour of the exact correlator and the truncated OPE for different numbers of power 
corrections retained in the expansion and for several values of the momentum transfer q. Clearly, retaining a fixed 
number of power corrections and requiring the accuracy of the truncated OPE to be better than say 1% for wT < 1.2 
restricts the region of the momentum transfer to "not very large" values. As is clear from Figs. [1] and [H in the region 
of uiT where the truncated OPE provides a good description of the exact correlator, the contribution of the excited 
states is still rather large, i.e., one is still rather far from the plateau. Therefore, a direct determination of the form 
factor from a truncated OPE is not possible. The procedures of the method of sum rules are aimed at modeling the 
contribution of higher states to the correlator and at obtaining in this way the ground-state form factor. 

5. SUM RULES AND EXACT EFFECTIVE CONTINUUM THRESHOLD 

The sum rules for the two-point correlator II and the three-point correlator T are merely an expression of equality 
of these correlators calculated in the "quark" and in the hadron basis: 

i? ge -^ T + n cxcitcd (T) = n (T) + n powor (T), (5.1) 

F 1 g (g)i? g e-^ T + r oxcitcd (T, (Z ) = T {T, q) + r powcr (T, q). (5.2) 
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Tope (T, g) r p E (T, q) r 0PE (T, q) 

r(T,q) r(T.q) T(T, q) 
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(a) (b) (c) 

Fig. 2: The accuracy of the truncated OPE for T(T, q) depending on the number of retained power corrections, for different 
values of the momentum transfer q. The plots present FopE(r, q)/T(T, q) for the following cases: (a) the terms including 0(T 6 ) 
are retained in Tope; (b) the terms including 0(T 8 ) are retained in Tope; (c) the terms including 0(T 10 ) are retained in Tope- 
Solid (red) line: q = 0; dashed (green) line: q = 1.5^/mtj; dash-dotted (blue) line: q = 3\Zmu). The vertical line at ujT = 1.2 is 
the upper boundary of the Borel window. 



Making use of the standard quark-hadron duality assumption — that the contribution of the ground state is dual to 
the low-energy region — we obtain the following duality relations: 1 

R g e- E * T = J dze" zT / 9o(2)+n powcr (T) = n dua i(T,z c n ff (T)), (5.3) 
o 

z ei! (T,q) z a «(T,q) 

F g (q)R g e- E z T = J dzi J dz 2 e~i zlT e-^ Z2T Ao(z 1 ,Z2,q)+T powcl .(T,q) = T duai (T,q,z c f f (T,q)). 



(5.4) 

We call the r.h.s. of (|5 . 3|) and (|5.4p the dual correlators Uduai(T, z^} s (T)) and Fduail!? 1 , q, z e s(T, q)), respectively. Let 
us notice that the dual correlators have both an explicit dependence on T and an implicit dependence on T via z e ff- 
As the consequence of (|4.4[) . n powcr (T) = r powor (T, q = 0). It is natural to set 

z^(T)=z eS (T,q = 0). (5.5) 

Then, due to the Ward identity (|4.14[) , the relations (|5.3[) and (|5.4[) yield the correct normalization of the form factor 
F g (q = 0) = 1, as required by current conservation. 

We treat the expressions (|5.3p and (|5.4[) as the definitions of the exact T- and g-dependent effective continuum 
thresholds z c s(T, q) and z^} s (T) corresponding to the true ground-state parameters on the l.h.s. The full T- and q- 
dependences of z e g (T, q) can be obtained by solving Eqs. (15. 3p and (|5.4|) for the known exact bound-state parameters 
i? g and F g (q) and the exact power expansions n power (T) and r power (T, q). In the HO model this can be easily done 
numerically. Without loss of generality we set m = ui and show the corresponding results in Fig. [3] Obviously, 
the effective continuum threshold does depend sizeably on both T and q. Now, a good or a bad extraction of the 
ground-state parameters depends on our capacity to find a reasonable approximation to the exact effective continuum 
threshold. 

Before closing this Section, we would like to make a comment on the T-dependence of the effective continuum 
threshold and the analytic properites of the dual correlator. 

Recall that a correlator in perturbation theory satisfies the standard dispersion representation in E with the cuts 
given by the quark diagrams. In the language of hadronic states, this correlator contains contributions of an infinite 
number of such states. 



The standard duality assumption may be formulated as follows: the contribution of higher hadron states is dual to the high-energy 
region of the free-quark diagrams. Therefore, excited states do not receive any contribution from power corrections. One may, of course, 
argue whether this assumption is physically relevant: it would be perhaps more natural to relate to the excited states also high-energy 
contribution of diagrams containing the confinement interaction, like those corresponding to T?2- This would, however, require a detailed 
knowledge of the confining potential. 
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Fig. 3: The effective continuum threshold z c g(T, q) for the 3-point function, obtained by solving numerically Eq. (|5.4[) using 
the exact bound-state parameters 7? g and F s (q) as well as the exact power expansion V povjel (T,q), versus the Euclidean time 
T at fixed values of the momentum transfer q (a) and versus q at fixed values of T (b). In (a) the vertical dashed lines identify 
the fiducial range in T (see text), while the solid lines are linear fits of z e ff(T, q) in the fiducial range. 



The dual correlator is a hand-made object and its analytic properties are very much different: the dual correlator 
should reproduce the contribution of a single ground state only, so it should have a single pole in E and no cuts. 
Equivalently, the Borel-transformed dual correlator should contain only single exponential exp(— E S T). Obviously, 
this may happen only if the effective threshold has a complicated dependence on T, since the spectral densities of the 
dual correlator and the perturbative correlator coincide. 

Therefore, the dependence of the effective threshold on T is its inherent property required by the analytic properties 
of the dual correlator. 



6. NUMERICAL ANALYSIS 



Let us now consider a restricted problem: assume that we know exactly the energy of the ground state, E g . Our 
goal is to determine the decay constant and the form factor of this state from the sum rules (|5.3|) and (|5.4[) . This 
problem is similar to the realistic case of the B-meson observables: in this case one knows the mass of the state and 
attempts to calculate its characteristics. 

To obtain theoretical estimates for ground-state parameters, we must perform the following steps: 

First, according to [l[ we should determine the Borel window (or the fiducial range), where according to the sum- 
rule philosophy the method may be applied to the extraction of the ground-state parameters: (i) the lower boundary 
of the T-window is found from the requirement that the ground state gives a sizable (we require more than 50%) 
contribution to the correlator; and (ii) the upper boundary of the T-window is obtained from the requirement that 
the truncated OPE gives a good (we require better than 1%) approximation to the exact correlator. 

Second, we must choose a criterion to fix our approximation for the effective continuum threshold z c g(T,q) and 
obtain in this way our sum-rule estimates for hadron observables. 

Third, we should define the way to obtain the band of values which contains the actual value with 100% probability 
and a flat probability distribution, in other words, define the way to obtain systematic errors of the hadron parameter. 
Obviously, this is a highly nontrivial point. 

The standard procedure adopted by sum-rule practitioners is to assume a T-independcnt effective continuum thresh- 
old. The quantity may be then either chosen as a q-independent constant or adjusted for any value of q, separately. 

The variation of the extracted form factor in the window is conventionally treated as the systematic error of the 
ground-state parameter. It should be recalled that in QCD this is a matter of belief: it is not possible to prove that 
the interval obtained in this way indeed contains the true value. 

In a series of recent publications we have demonstrated that in many cases these standard criteria fail: the actual 
parameters of the ground state turn out to lie beyond the range predicted by the sum-rule procedures. In [ij we have 
shown this for the extraction of the ground-state decay constant from the two-point function; in [5J, [6| we have given 
examples of the form factor at specific values of momentum transfers, for which the standard criteria (stability in 
the window, deviation of the energy of the cut correlator from the ground-state energy in the window) suggest an 
accurate extraction of the form factor, whereas in practice the actual error turns out to be much larger. Now, we 
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are going to investigate whether improvements may be obtained by releasing the requirement of the T-independcnt 
effective continuum threshold and by imposing other criteria than Borel stability for obtaining the error estimates. 
In this work we discuss three different Ansatze for the effective continuum threshold: 

z eS (T,q) « z c (q), (6.1) 
Z eS (T,q) « z^(q) + z^(q) uT, (6.2) 
z cff (T,q) » z$(q) +z?(q) toT + z® (q) co 2 T 2 . (6.3) 

At each value of q we fix the unknown parameters on the r.h.s of Eqs. (|6.I[) - (|6.3[) in the following way: we define the 
dual energy, E dual (T,q), as 

EduoL(T,q) = --^ log T dua i(T,q,z cS (T,q)), (6.4) 

where r^aiC? 1 , q, z D g) is the r.h.s. of Eq. (|5.4|1 calculated using the approximations (|6.ip - ()6.3p for z e ff(T, q). Let us 
emphasize that the implicit dependence on T via z e g is crucial for the calculation of the energy: if the exact effective 
continuum threshold has a sizeable dependence on T, neglecting this dependence leads to a completely wrong energy, 
which we later tune to the exact known value. This is, in fact, the main source of error in the sum-rule predictions 
for hadron observables. 

Next, we calculate Ed na l(T, q) at several values of T = Tj (i = 1, . . . , N) chosen uniformly in the fiducial range. 
Finally, we minimize the squared difference between the dual energy Eduai(T, q) and the exact value E g : 

N 

,2 



X ~ N 

i=l 



l J J2[E dn&1 (T i ,q)-E s } 2 . (6.5) 



A. Decay constant 



Let us start with i? g , the square of the decay constant. The standard procedures for this case were discussed in 
detail in The main results of Q may be summarized as follows: in spite of the fact that the numerical value of the 
extracted decay constant turns out to be not far from the exact value (see also the discussion in [14j). the standard 
procedures fail to produce realistic error estimates: namely, the known true value lies outside the range obtained from 
the standard sum-rule analysis. 

Fig. |4] shows the results obtained by allowing for a T-dependent threshold and tuning its parameters according to 
(|6.5p : one can see three approximations for z^ s (T) obtained by minimizing the function (|6.5[) for 

£dual(T) = ~ log II dual (T, z c n ff (T)), (6.6) 

riduai(r, zj^) being the r.h.s. of (|5.3p . and the corresponding Rdu&\(T)/R g . The extraction of R is presented for three 
cases where a different number of power corrections — three, four, and infinite (i.e., exact power corrections) — in 
IIope are taken into account. 

The case of the exact power corrections is presented as an illustration: we extract the ground-state parameters from 
the correlator by applying the sum-rule machinery in the window 0.7 < wT < 1.2 (although in this case, of course, 
the parameters may be extracted directly from the large-T behaviour of this correlator) . 

In all cases one can see an obvious improvement when the linear Ansatz z£ s (T) instead of the constant z^ s is used. 
Going beyond the linear approximation and allowing for the quadratic function z® s (T), however, does not lead to 
further improvement in the realistic cases of three and four power corrections in IIope- This effect is just the reflection 
of the fact that the ground-state parameter may be extracted from the correlator only with a limited accuracy, even 
if the ground-state mass is known precisely (see also the detailed discussion in [4|). The same situation occurs for the 
form factor in Section [6 Bl 



B. Form factor 

The standard sum-rule analysis of the ground-state form factor was presented in || . We have given there an explicit 
example of the sum-rule extraction of the form factor at a specific value of the momentum transfer where the method 
clearly fails in the following sense: the true value of the form factor turned out to lie well beyond the interval obtained 



9 



z=ff(T)/6j E du „(T)/E g Rd„ai(T)/R g 




0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 0.6 0.8 1 1.2 0.6 0.8 1 1.2 



Fig. 4: The extraction of R s , the square of the decay constant of the ground state, making use of TV power corrections in 
IIope(T): N = 3 (first row), TV = 4 (second row), and N = oo (third row). First column: the constant, linear, and quadratic 
approximations for the exact effective threshold z^} s (T) obtained by minimizing (|6.5I) . Second column: the corresponding 
-Eduai(T) / E e . Third column: the corresponding fldual(T) / R g . Red (dotted): constant z^g', blue (dashed): linear z^ s (T); green 
(solid): quadratic z® B (T). The vertical lines indicate the boundaries of the window. 

by the standard sum-rule procedures. Moreover, the actual errors in the case of the form factor were shown to be 
much larger than those for the case of the decay constant. 

We now apply a T-dependent Ansatz for the effective continuum threshold to extract the form factor. Recall that 
because of current conservation, the form factor F(q) should obey the absolute normalization F(q = 0) = 1. We 
therefore require z^ s (T) = z e g(T, q = 0): then at q = the r.h.s. of Eq. (|5.4|) reproduces the r.h.s. of Eq. (|5 ,3[) and 
the form factor is automatically properly normalized. Therefore, the form factor extracted from the sum rule at small 
momentum transfers is anyway not far from its true value. Of real interest is the extraction of the form factor at the 
intermediate momentum transfers q/uj = 1 -j- 2. 

In the HO model we know all power corrections exactly, therefore, we have no limitation on the upper boundary 
of the window. However, to be in line with the realistic situation where only a limited number of power corrections 
is available, we define the window to be 0.7 < ioT < 1.2 (see Fig. [2]). 

The results for the form factor F(q) obtained from the dual correlator by optimizing the parameters of the three 
Ansatze (|6.1[) - (16.3[) according to (|6.5[) are shown in Figs. [5]and[6l For a T-independent approximation (|6.1|) , the 
minimization of \ 2 (|6.5|) leads to Zq (q) for which the dual energy E^ ua y(T, q) differs from E g by less than 1% and 
the dual form factor i"d U ai(<7) turns out to be practically T-independent in the whole fiducial range. Such a stability, 
usually referred to as the Borel stability, is often (erroneously) claimed to be the way to control the accuracy of the 
extracted form factor. From Fig. [5] it is clear, however, that the actual error of the extracted form factor turns out to 
be much larger than the variation of the form factor in the fiducial range of T. Moreover, the T-independent Ansatz 
is applicable only for not very large q: for q > 1.7 u) it fails to reproduce the ground-state energy in the fiducial range. 

Allowing for a linear Ansatz (|6.2p for z e g (T, q) extends the range of the momentum transfers where the form factor 
can be extracted from the sum rule with a reasonable accuracy. Note that the exact effective continuum threshold 
of the HO model can be well approximated by a linear function of T in the whole fiducial range, see Fig. [3£a). 
Nevertheless, this does not mean that the minimization of (|6.5p finds the exact effective continuum threshold and 
leads to the exact form factor: the deviation of the extracted form factor from the exact one amounts to a few percent 
for q < 2u> and increases dramatically for q > 2 u>, see Fig. [51 The dramatic increase of the error at q > 2u is directly 
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Fig. 5: The extraction of the ground-state form factor at q/uj — 0.5 (first row), q/uj = 1.5 (second row), and q/u) = 2 (third 
row). First column: the effective continuum threshold as found by the fitting procedure (|6.5[) ; second column: the fitted dual 
energy Eduai(T,q)/E s ; third column: the corresponding form factor ratio Fduai(T,q)/ F g (q). Dashed black line: exact z e s', full 
circles (red): constant z^ s \ empty squares (blue): linear z^ s (T); full diamonds (green): quadratic zf s (T). 
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related to the fact that the contribution of the ground state to the correlator decreases rapidly with q in the given 
T-window 0.7 < 10T < 1.2 (see Fig. [T]). Let us emphasize that the actual error of the extracted form factor turns out 
to be much larger than the variation of the form factor in the window, see Fig. [5] 

One may try to go further and to consider the quadratic Ansatz (|6.3[) . From Fig. [5] it is clear that the quadratic 
Ansatz leads to certain instabilities in the extracted value of the form factor and does not, in general, improve the 
actual accuracy of the form factor extraction. This is not strange: such a behaviour just reflects the fact that the 
sum-rule extraction of the ground-state parameters can be performed only with a limited accuracy, which has been 
known since the initial SVZ paper [l| . It should be therefore clear that there is no way to get further improvement in 
the accuracy of the extracted form factor by increasing the degree of the polynomial approximation of z e g(T, q). 

Nevertheless, we have seen a clear improvement in the outcome of the extraction procedure as one goes beyond 
the assumption of a T-independent effective continuum threshold: first, the actual accuracy turns out to be (much) 
better; second, the T-dependent Ansatz allows one to extract the form factor in a broader range of the momentum 
transfer. 

The crucial question for the application of sum rules to the bound-state form factors is how to understand the 
actual accuracy of the ground-state parameter extracted from the sum rule. 

(i) Obviously, the Borel stability does not guarantee the good extraction of the ground-state parameter and the 
variation of the extracted ground-state parameter in the window does not provide a realistic error estimate. 

(ii) It should be understood that even a very accurate reproduction of the ground-state energy in the window does not 
lead to equivalently accurate extraction of the ground-state parameters. This is rather unpleasant since the deviation 
of the dual energy from the known ground-state energy is the only visible characteristic which can be controlled in 
the realistic cases. What really matters for the extraction of the bound-state parameters is the deviation of the fitted 
effective threshold from the exact effective threshold which, however, in the realistic cases remains unknown. 

In the HO model, a realistic error estimate of the accuracy in a wide range of values of q may be obtained by 
comparing to each other the form factor values extracted by assuming the linear and quadratic Ansatze for the 
effective continuum threshold (|6.2p and (|6.3|) . Whether this feature persists in QCD is an interesting and important 
issue to be addressed in the future. 



7. CONCLUSIONS 

We studied the extraction of the ground-state form factor from the vacuum-to-vacuum correlator in the exactly 
solvable harmonic-oscillator model applying the standard procedures of the method of QCD sum rules. Let us 
summarize the main messages of our analysis: 

• The knowledge of the correlator in a limited range of relatively small Euclidean times T (that is, large Borel 
masses) is not sufficient for the determination of the ground-state parameters. In addition to the OPE for the 
relevant correlator, one needs an independent criterion for fixing the effective continuum threshold. 

• Assuming a T-independent (i.e., a Borcl-parameter independent) effective continuum threshold, the error of the 
extracted form factor F(T, q) turns out to be typically much larger than (i) the error of the description of the 
exact correlator T by the truncated OPE Tope and (ii) the variation of F(T,q) in the fiducial range of T (i.e., 
the Borel window). As the result, the actual value of the hadron form factor may lie far outside the range 
covered by the sum-rule estimate F(T, q) within the Borel window of T. Therefore, the variation of the hadron 
form factor (as well as other hadron parameters) within the Borel window cannot be used as an estimate of 
its systematic error. The latter point is of particular relevance for the practical applications of sum rules in 
QCD, since the Borel stability is usually (erroneously) believed to control the accuracy and the reliability of the 
extracted ground-state parameter . 

• In the cases where the ground-state mass is known (e.g., experimentally measured), the actual accuracy of the 
sum-rule analysis may be considerably improved by allowing for a T-dependent effective continuum threshold 
and finding its parameters by minimizing the deviation of the energy of the dual correlator from the known 
value of the ground-state energy, Eq. (|6.5[) . 

In the HO model, this procedure was shown to yield clear improvements in the extracted values of the decay 
constant and the form factor. 

Moreover, in the HO model the deviation between the hadron parameter values obtained by using different 
Ansatze for z e s gives de facto a realistic error estimate. This was observed for both the decay constant and the 
form factor. 



12 



Unfortunately, it remains impossible to construct the band of values which may be proven to contain the actual 
form factor. In this sense, the method of sum rules is not able to provide rigorous error estimates. Nevertheless, the 
application of the new procedures formulated in this paper to hadron form factors in QCD seems very promising with 
respect to improving the actual accuracy of the method. 
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